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We study a topological updating scheme in three dimensional U(l) gauge theory. Some expectations for four 
dimensional SU(N) gauge theories are discussed. 



1. Introduction 

The presence of topological objects in a 
lattice theory may play havoc with updating 
schemes. The problem is that changes in the 
number of these objects usually involve fairly 
large changes in quite a number of lattice vari- 
ables, resulting in large changes in the action, 
and hence in a very low acceptance rate. In 
other words, there is a very poor sampling over 
the different topological sectors, characterized by 
the number of said objects. Not only in theo- 
ries that contain objects of a genuine topolog- 
ical nature, such as instantons and topological 
charge, do these problems arise, but also in the- 
ories that possess quasi-topological objects, such 
as monopoles and vortices. 

A strategy to tackle this problem is to devise 
an algorithm that proposes these changes with a 
large probability of acceptance. The cluster al- 
gorithm succeeds in doing just that for, e.g., the 
XY models[l]. When this proves to be too dif- 
ficult, e.g. for SU(N) gauge theories in four di- 
mensions, we may consider 'adding' transitions 
by hand[2]. The question here is how to get a 
reasonable acceptance probability for these 'topo- 
logical updates'. By reasonable we mean larger 
than, say, 5%. We will address this problem in 
3-d U(l) pure gauge theory. 

2. The model and the method 

The Wilson action for 3-d U(l) gauge theory 
is given by 

S = (3 J2 ReU^A^) (1) 



In this theory, we can define monopoles in the 
following way. Writing Ufi^(x) = exp(i 61^^ (a;)), 
we can split 6 : 9fj,v{x) = 9 fj,v{x) -\- 27rn^y(a;), 9 G 
< -TT, tt], n integer. Note that n = 0,±1, ±2. We 
can transpose this to the dual lattice by — 
/i) = pXii, p{x) , which satisfies the following 
relation (dropping the hat from now on): 

m{x) = 0,±1,±2,±3. (2) 

When m(x) > 0(< 0), there is a (anti-) monopole 
at site X on the dual lattice. When m(x) = 0, 
but any of the six j's connected to x is non-zero, 
a Dirac string passes through x. It is easy to see 
that then an even number of j's are non-zero. 
Nothing prevents j from forming closed loops. 
This is indeed what is found [3]. We can now dis- 
tinguish between two types of loops: ones that 
are closed through the boundary, and ones that 
are not. Although we can assign a winding num- 
ber per direction to former type, it turns out that, 
in each direction, only the total sum of winding 
numbers is gauge invariant. It is possible to split a 
single loop with winding number 2 into two loops 
with winding number 1, merely by performing a 
gauge transformation. 

These Dirac strings introduce meta-stability[4, 
5]. Suppose we have a gauge configuration with 
a single Dirac string with winding number 1. In 
order to change this number, we need to create 
another string. As this involves changes over large 
distances, the probability of this happening will 
be (very) low in, e.g., a standard Metropolis al- 
gorithm. 

If we now take a slice of the lattice through 



which a Dirac string pierces, and look upon this 
as a two dimensional gauge configuration, we find 
that the topological charge of this 2-d configura- 
tion is equal to the winding number of this string, 
i.e. ±1. So we now have a scheme to make a 
topological update for 3-d U(l): we simply add a 
stack of 2-d U(l) configurations with topological 
charge (or -1) to our 3-d U(l) configuration, 
and perform a global Metropolis on this new con- 
figuration. 

There are different types of 2-d configurations 
we can use. We can use configurations from an 
equilibrium ensemble at a given /32-d- '-^^ 
also construct gauge configurations with constant 
field strength. We will use 2-d gauge configura- 
tions from several /32-d values, as well as specially 
constructed configurations: 



Ui{x) = exp(iw (a;2 - 1)), 
U'jix) = 1 for X'j = I, ■ ■ 

Ui{x) = exp(iw (a;2 - 1)), 
U2(x) = exTp(—iuiNxi) 
for X2 = N 



,{N-l) 



(3) 



where oj = 2i: / N'^ . This leads to a gauge configu- 
ration with a field strength of Up(x) = exp(— iw) 
throughout the lattice. The topological charge 
for this configuration is Q = —1. 

From a theoretical point of view we can make 
an estimate of the acceptance rate of these topo- 
logical updates. From [6], we find 

n = erfc{^y^), hi=<AS>, 
where AS = S(U^,)- S(U^, + U2.d) . (4) 
From this, we expect the following. 



type 


AS 


n 


eq. (3) 
2-d conf. 


/3#3 <:Up> [1- < >] 


0.252 



for a 16^ lattice at /3 = 2.50, for which 
< Up >=0. 8542(6). < U^,'' > >0.9 for P^'' > 4.0. 

3. Results 

So from this analysis we expect no topologi- 
cal updates for any 2-d gauge configurations but 



/? 


acc. rate 


<Up> 


Q 


M 


1.50 


0.436(42) 


0.6880(5) 


0.4249 


222(2) 


1.75 


0.374(21) 


0.7619(4) 


0.3644 


75.0(1.5) 


2.00 


0.334(37) 


0.8061(3) 


0.3186 


22.7(1.0) 


2.25 


0.316(21) 


0.8342(5) 


0.2819 


6.6(1.0) 


2.375 


0.164(44) 


0.8453(2) 


0.2658 


1.9(2) 


2.50 


0.152(45) 


0.8542(1) 


0.2516 


0.97(15) 


2.75 


0.128(45) 


0.8700(1) 


0.2244 


0.21(2) 


3.00 


0.096(33) 


0.8809(1) 


0.2017 


0.054(8) 


3.25 


0.032(12) 


0.8909(1) 


0.1814 


0.015(4) 


3.50 


0.056(31) 


0.8992(1) 


0.1635 


0.004(2) 


3.75 


0.040(25) 


0.9036(1) 


0.1482 




4.00 


0.024(16) 


0.9125(1) 


0.1335 




4.25 


0.024(10) 


0.9180(1) 


0.1208 




4.50 


0.020(8) 


0.9228(1) 


0.1095 




4.75 


0.012(16) 


0.9270(1) 


0.0994 




5.00 


0.008(11) 


0.9309(1) 


0.0902 





Table 1 
value. 



The results for 5k measurements per /3 



the one with a constant field strength. This is 
indeed what is found. We made 15k updates, be- 
ing 5-hit Metropolis, after discarding 2k updates. 
Dirac strings were proposed every 10 updates. We 
indeed found that no topological updates were ac- 
cepted for any 2-d gauge configurations but the 
one with a constant field strength. Here, we found 
an acceptance rate of 0.133(17), which is remark- 
ably smaller than the 'expected' 0.252. 

Returning to 2-d gauge configurations from an 
equilibrium ensemble; we might consider choos- 
ing a very large /3'^"*^ . This will lead to an average 
plaquette, which will be very close to unity, hence 
rendering a reasonable acceptance rate. However, 
since such a two dimensional gauge field is re- 
quired to carry a topological charge Q = ±1, we 
will find that all plaquettes but one or two will 
be very close to unity, while the remaining one or 
two will be close to minus unity. For these con- 
figurations we get AS = 4/3 < Up >= 8.542 and 
thus < 10-8. 

Being surprised by the large deviation of the 
acceptance rate we found at /3 = 2.50 from the 
theoretical one, we made shorter runs at several 
/3 values to check the validity of our formula (4). 
The results are as follows: 



The transition around [3 = 2.3 is striking. It 
would appear that Dirac strings become sup- 
pressed for larger [3 values. This becomes 
very reasonable when we consider the average 
monopole number, M =< |m(a;)| >, pre- 
sented in the last column of table 1. 

So the transition in the acceptance rate occurs 
where the average number of monopoles in a con- 
figuration drops below 2. It is then not at all 
surprising that Dirac strings also become (very) 
rare. 

4. Outlook for SU(N) 

The requirement for topological updates to 
have a reasonable acceptance rate turns out to 
be a very restrictive one on the actual implemen- 
tation of the algorithm. As a matter of fact, only 
the smoothest updates possible stand a chance 
of having a reasonable acceptance rate. We also 
found that there can be large deviations from 
the theoretical predictions, due to the intrinsic 
properties of the model, which have not been ac- 
counted for in the analysis. 

Let us now speculate on 4-d SU(N) gauge the- 
ories. As all these theories allow for topological 
charge and instantons, we expect very large tun- 
nelling times from one topological charge sector 
to another for sufficiently large gauge coupling [3. 
As we have seen, taking a realistic gauge configu- 
ration will lead us nowhere, so we have to resort 
to specially constructed gauge configurations. 

Taking SU(2) for the moment, we can construct 
two types of configurations. The first one has 
constant field strength for each {jJ.ji') plane, and 
is given by 

U^{x) = exp(ir 

-LyUj^,yX^,8^^^(L^.i)]) (5) 

where, from the co-cycle condition = 
27rn^y/_L2, n^^ integer. Note that n^^ 7^ n^^ 
is allowed, r is one of the Pauli matrices, or 
possibly a properly normalized linear combina- 
tion of them. This configuration carries topolog- 
ical charge Q = 2[(ni2 - n2i)(n34 - n43) - (nis - 



n3i)(n24 - n42) -|- (ni4 - n4i)(n23 - n32)] Notice 
that Q is even. The other configuration is the 
well-known single instanton configuration [7-9]. 

These configurations all have rather large ac- 
tions. The one with constant field strengths has 
an action of S/(3 = 2-17^ Y.^^ (n^^ - ny^)^, while 
the single instanton configuration has an action of 
S/ [3 = 9.6, but which can be lowered to S/ [3 = 6.8 
by cooling[8]. 

These numbers are not encouraging from the 
topological update's point of view, but it should 
be kept in mind that for non-Abelian theories, 
< Ai5 > may turn out to be smaller than ex- 
pected from this analysis. 
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